Experimental investigation of optical atom traps with a frequency 

jump 

P Ahmadi, G Behinaein, B P Timmons and G S Summy 
Department of Physics, Oklahoma State University, Stillwater, Oklahoma 74078-3072 

Abstract 

We study the evolution of a trapped atomic cloud subject to a trapping frequency jump for two 
cases: stationary and moving center of mass. In the first case, the frequency jump initiates oscilla- 
tions in the cloud's momentum and size. At certain times we find the temperature is significantly 
reduced. When the oscillation amplitude becomes large enough, local density increases induced by 
the anharmonicity of the trapping potential are observed. In the second case, the oscillations are 
coupled to the center of mass motion through the anharmonicity of the potential. This induces 
oscillations with even larger amplitudes, enhancing the temperature reduction effects and leading 
to nonisotropic expansion rates while expanding freely. 

PACS numbers: 32.80.Lg, 32.80.Pj 
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I. INTRODUCTION 



A trapped atomic cloud subjected to a sudden change in its trapping potential can ex- 
hibit a wide range of behavior with many possible applications. It can be used as a model 
for studding fundamental phenomenon ranging from superfluidity to the generation of non- 
classical states of a matter- wavepacket. For instance, a Bose-Einstein condensate can be 
reversibly created by suddenly changing the trapping potential sinusoidally moving a 
magnetic trap center can be used to observe the excitation modes of a Bose-Einstein con- 
densate Q; the modulation of optical traps can be used for coherent control of the center 
of mass wavepacket ;3|; and the breathing mode oscillation of atomic clouds in ID and 3D 
optical lattices can be real i Z e d by sudden ehan ging the la «ce depth Q. The later has a, SO 
captured a great deal of theoretical interest because of the possibility of creating squeezed 



states of a matter wavepacket 



m 



ul- 



na. 

The oscillation of an atomic cloud after a sudden jump in the trapping frequency has also 
been proposed as a tool for optical cooling and hence controlling the onset of Bose-Einstein 
condensation (12! Q] . This cooling method can be categorized among the techniques that 
use geometrical manipulation of an external potential to reduce the temperature of atoms 
or molecules. For example, Flories et al. jl4] used a spatially varying electric field to 
longitudinally cool a molecular beam to 250 nK, while Ketterle and coworkers achieved pico- 
Kelvin temperatures by adiabatically decompressing a gravito- magnetic trap ^5|. Exposing 
a freely expanding cloud of cold atoms to a pulsed potential has also been considered as a 
method for temperature reduction |l6j ]. 

In this paper we use optical traps to experimentally realize a potential with a sudden 
frequency jump so as to explore the cooling proposals of Ref. [l^. These traps have 
become versatile tools for atom optics research, with the observation of all-optical Bose 
Einstein condensation an example of their dramatic potential [l^. Early ex per iments on 
intensity modulation of the laser beam used to form an optical trap JJ, lisl. IliJ. I20I ] were 
used to demonstrate the parametric excitation of an atomic cloud. In this paper we will 
discuss experiments with these traps, comparing the observations with theory, and presenting 
possible future applications. The format of this paper is as follows. In Section II the 
experimental configuration is reviewed. In Section III we first lay out the theory of non- 
adiabatic cooling of an atomic cloud, then follow by presenting experimental data for the 
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case of a sudden frequency jump with a stationary center of mass. This discussion is then 
repeated in Section IV for the case of a moving center of mass. 

II. EXPERIMENTAL CONFIGURATION 

Our experimental apparatus has been described previously (21], |^ so only a brief de- 
scription is given here. We create a potential for Rb87 atoms using their interaction with 
10.6 fim light from a CO2 laser. The CO2 laser radiation frequency is far enough below 
that of the resonances of the atoms that its electric field can be considered quasi-static. 
This induces a dipole moment proportional to the local electric field and a lowering of the 
atomic ground state energy given by, U = —^a g \E\ 2 , where a g is the ground state static 
polarizability and \E\ 2 is the time averaged square of the laser light's electric field. 

Two CO2 beams were directed into a vacuum chamber in a geometry such that they 
propagated orthogonally to each other. One of the beams propagated in the vertical direction 
(x direction) and the other in the horizontal direction (z direction). Crucially the position 
where the x beam intersected the z beam could be adjusted. That is, the foci of the two 
beams did not necessarily coincide. The light for these beams originated from a 50 Watt RF 
excited CO2 laser whose total power was controlled by passing the output light through an 
acousto-optic modulator (AOM). The first order beam of the modulator was then directed 
into another AOM which was used as a beamsplitter whose ratio could be changed by varying 
the input RF power. This ratio was used to control the time dependent potential. In order 
to fulfill the conditions necessary for this experiment, the CO2 beams were aligned such that 
the x beam was the first order of the second AOM and the z beam the zeroth order of the 
same AOM. The data were taken by destructively imaging the cloud using a resonant probe 
laser which passed through the atom cloud and was then incident on a CCD camera. A 
gaussian fit to the optical density data determined the spatial extent of the atomic cloud. 

III. STATIONARY CENTER OF MASS 

We begin by theoretically considering the situation where the center of mass of the atomic 
cloud remains stationary after a sudden change in the trapping potential. We take a sample 
of atoms in thermal equilibrium inside a harmonic potential of frequency u)q which is non- 
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adiabatically changed to oj\ at t = 0. Neglecting the effect of collisions and following Balatov 



et al. 



12j, the momentum and position probability distribution widths are given by, 



"lit) = o^(0)[l + i-f + (1 - (-) 2 ) 008(^)1 (1) 



k(o)[i+(^) 2 +(i-(^ 

0*(t) = k(0)[l + (^) 2 + (1 - (^) 2 ) 008(^)1 (2) 



where at t = ^- they reduce to, 



= (^)VJ(O). (4) 

For c^! <C ujq, a narrow momentum distribution is produced corresponding to a lower effective 
temperature. In the case of a single Gaussian laser beam propagating in the z direction the 
average electric field can be expressed as, 

\E(x,y, z)\ 2 = ^^exp ~ 2(x2 +/ 2) (5) 
w(z) z w(z) z 

where w(z) = wq (l + (^) 2 ) 2 an d Eo is the electric field amplitude. Here w is the beam 
waist at the focus, and is the Rayleigh length. By carrying out a series expansion around 
x = y = z = and discarding terms of third order and higher, a harmonic approximation 
of the potential is obtained, 

U(x,y,z) = U ( 2{x2+ 2 y2) +^) (6) 
V w o z r) 



where U = —^a g EQ. The oscillation frequencies in the three directions are u% = u 2 ~ 
AUo/mvjQ,uj 2 z ~ 2Uo/mz R . When the total laser power is shared with a second beam propa- 
gating in the x direction which crosses the first beam at its focus, these frequencies change 
to ujy ~ 4Uo/mwQ, ijj 2 x = u>1 ~ 2Uo/mw 2 l at the intersection point. This assumes that the 
beams have identical properties and zr 3> Wq. According to these results, an abrupt change 
from a two beam to a one beam geometry will produce a significant frequency change along 
the z direction. The ratio of one and two beam frequencies is, 

^lbeam _ ^0 _ ^ ^ 
W 2 beam Z R 

which is independent of the total laser power and is determined by the beam waist at 
the focus. Since the temperature of the atomic cloud is proportional to the momentum 



distribution squared, combining Eqs. © and (J7J) gives 

Tfinal _ /^A 2 _ ( A 
initial \U Q ) 

To explore this behavior experimentally we aligned the CO2 beams such that their foci 
overlapped. Atoms were loaded into the potential minimum created by the intensity maxi- 
mum located at the intersection of the beams. The power in the x beam was transferred to 
the z beam by switching the RF power on the second AOM to zero. This abrupt transfer of 
power mostly produced a change in the effective frequency in the z direction and initiated 
the oscillation of the trapped cloud. If the shift of power to the z beam happens immediately 
after loading, the FORT will not produce the desired oscillation since the atomic density is 
approximately 10 14 atoms cm -3 and cold binary scattering dramatically affects the evolution 
of the cloud. Therefore the collisionless approximation of the previous section can not be 
used. To move into the collisionless regime it is necessary to decrease the atomic density by 
reducing the number of atoms before the power transfer. This is accomplished by reducing 
the total power in the CO2 beams in order to carry out forced evaporation of hot atoms. 
This reduces both the density of atoms in the trap and their temperature. This process is 
carried out so that the atomic cloud which remains in the FORT has densities of around 10 11 
atoms cm -3 . This guarantees that the mean free path of the individual particles is bigger 
than the cloud size so that the evolution of the sample can be considered in the collisionless 
regime. 

After evaporative cooling the total power was transferred to the z beam and the cloud size 
started to oscillate in the z direction. However, the time evolution of the cloud size did not 
obey the simple periodic form predicted by Eq. (j2J. Our observations show that a domain 
with high atomic density appears at the center of the trap periodically during the cloud's 
evolution . This structure is reminiscent of the two component structure which appears at 
the onset of a degenerate Bose gas. Three cross sections of the cloud along the z-axis are 
shown in FigO Figur^(a), (b) and (c) were taken 2, 5 and 10 ms after switching off the x 
beam respectively. The total power in the CO2 beams was 5 Watts. A high density region 
at the center of Fig. E^b) is clearly visible. It disappears later in the evolution as seen in 
Fig. W^c). Creation of these high density regions is a direct result of the strong anharmonic 
potential and is predicted in jl3|. A Monte Carlo simulation of the anharmonic system can 
reproduce this behavior for the cloud. Figure 121 presents the result of such a collisionless 
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FIG. 1: The accumulated intensity along the function of z distance. Two beams from a 

CO2 laser crossed each other at right angles with overlapping foci and 5 Watts power, (a), (b) and 
(c) were obtained 2, 5 and 10 ms after extinguishing the x beam. A denser component is obvious 
at the center of the cloud in (b). 

calculation at the moment where the two component structure appears. The simulation was 
carried out for 3 x 10 4 atoms initially at thermal equilibrium with a 25/zm initial cloud size. 
It is worth noting that the high density regions never appeared when harmonic potentials 
were simulated. 

The high density region makes the observed optical density profiles complex enough that 
it is impossible to fit a Gaussian function to the data and infer the width of the cloud. 
Therefore to observe the oscillation of the cloud size, predicted by Eq. (J2J) we have reduced 
the total power in the CO2 beams to 0.5 Watts. In this regime the cloud temperature is low 
enough (below 2 /xK) that the oscillation amplitude does not move the atoms away from 
the region where the harmonic approximation of the potential is valid. Furthermore, as 
mentioned above, the evaporative cooling also reduces the atomic density to « 10 11 atom 
cm -3 , suppressing the two-body collision rate. The time dependency of the cloud size for 
this condition is given in Fig. El This shows that the cloud does not collapse to its original 
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FIG. 2: (a) Calculated z direction phase space distribution of the atoms at the moment where 
the two component structure appears. The atoms where initially distributed uniformly throughout 
the phase space and physically located at the cross over of two beams with overlapping foci. The 
horizontal axis is the z-direction spread of the particles in units of the beam waist and the vertical 
axis is the z-momentum in units of the photon recoil from a 780 nm photon, (b) The corresponding 
histogram of the particles' distribution in the z direction. 

extent during its evolution in the trap. We expect that this is the result of a small number 
of collisions. According to this data the size of the cloud in the z direction expands to about 
6 times its initial value. From Eq. (J3J) this implies a factor of 6 reduction in the width 
of the momentum distribution or a factor of 36 in the effective z-direction temperature. 
This number in combination with Eq. (JHJ) determines the beam waist at ~ 20/im, in good 
agreement with the theoretically calculated value of ~ 23/mi. 

IV. MOVING CENTER OF MASS 

In contrast to most other work with FORTs, we also studied geometries where the center 
of mass moves as a result of a change in the potential. This was achieved by offsetting the 
foci of the CO2 beams from each other. In this configuration, when one of the beams was 
switched off, the atoms moved towards the new potential minimum located at the focus 
of the remaining beam. For a harmonic trap the center of mass and internal dynamics of 
the atomic cloud are decoupled from one another. However, for an anharmonic trap the 
periodic motion of the center of mass of the atoms couples to the internal dynamics. Figure 
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FIG. 3: Evolution of the z cloud size inside the z beam for a stationary center of mass. Two beams 
from a CO2 laser crossed each other at right angles with overlapped foci. The power in the x beam 
was abruptly transferred to the z beam at t=0 to initiate the cloud's oscillation. 

@] schematically displays the trapping potential before and after power is transferred to the 
z beam. 

The observed z cloud size as a function of time for one period is presented in Fig. |S| As 
occurred in the previous section, the cloud size in this direction never shrinks to its initial 
value. Figure El also shows that the cloud undergoes a considerable increase in width by the 
time it reaches the focus of the z-beam at t ~ 10ms. It can be seen that the cloud size 
increases by as much as 10 times its initial value. This implies a compression in momentum 
spread by a similar factor and thus we expect a two order of magnitude reduction of the 
effective temperature in the ^-direction. This is significantly more than was possible in the 
stationary center of mass experiments. 

We also performed experiments in which the potential was switched off completely after 
the atoms had evolved for various times inside the z-beam. There are several features of 
the cloud's evolution after release from the z-beam that are worthy of note. Firstly, it is 
possible to tune the center of mass velocity in the ^-direction with high precision. This 
can be accomplished by changing the time span that the cloud evolves inside the z-beam 
before it is switched off. Such a property could be useful for q uantum reflection experiments 
where precise control of the impact velocity is required [23]. A second aspect of the free 
evolution is the observation of a collapse of the cloud in the ^-direction during the first 
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FIG. 4: Schematic of the dipole potential (a) before and (b) after transferring the power of the x 
beam to the z beam. In (b) the center of mass of the atomic cloud oscillates around the Gaussian 
focus of the remaining light, (c) cloud size after evolving inside the z beam for approximately half 
a period. The difference in velocities at the two ends of the cloud gives rise to a collapse after the 
potential is switched off for a TOF experiment 
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FIG. 5: The cloud size (closed circles) and center of mass position (open squares) in the z direction 
while evolving inside the z beam. Initially two beams from a CO2 crossed each other at right angles 
with their foci displaced by ~ 140 /mi. The power in the x beam was abruptly transferred to the 
z beam at t=0. 
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FIG. 6: Experimental data showing the cloud size. The power was 5 Watts per CO2 beam. The 
circular and square symbols represent the cloud size in the x and z directions respectively, (a) 
shows the results after various expansion times for the case where the FORT has been directly 
released from the crossed beams. Note that the circle and square symbols overlap each other in 
this case. In (b) the x and z cloud sizes are shown for the condition where the power from the x 
beam was switched to the z beam 15 ms before releasing the FORT. In (c) the ratio of the x to z 
cloud sizes of (b) is presented. 

few milliseconds after the atoms are released. To observe this collapse the cloud is first 
allowed to evolve inside the z-beam for nearly half a period of the center of mass motion. 
The collapse can be explained by the velocity difference between the atoms on either side 
of the cloud imprinted by the gradient of the trapping potential. This can be seen in Fig. 
IHc). The exact time at which to switch the z-beam off to observe the maximum collapse 
is determined experimentally. This time is 16 ms for the experiment shown in Fig. The 
resulting collapse or focusing can be used to transfer energy from the z to the x and y 
directions if the cloud has high enough densities such that the free expansion starts when 
the cloud is in the collisional regime. Under this condition, since at the release time the trap 
has become elongated in the z direction, the atoms moving along the z direction experience 
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a larger number of collisions compared to the x and y directions and part of their z kinetic 
energy transfers into x and y while expanding |24j . To explore this possibility, unlike in the 
previous section, the evaporative cooling is carried out such that the cloud had a density 
of ~ 10 14 atoms cm -3 (this corresponds to a mean free path of 9/im for our system). The 
velocity distribution should remain constant with time after the expansion makes the mean 
free path larger than the cloud size and the atoms reach the collisionless regime. Therefore 
;o determine the temperature the asymptotic expansion rate of the cloud data must be used 



25| after it has reached the collisionless regime . 



For comparison we have conducted another series of experiments in which the power in 
the CO2 laser was switched off without transferring total power to the z beam. Figure 
Efa) shows the size of the atomic cloud as a function of time after release from a crossed 
beam trap with 5 Watts of power. As can be seen, the expansion of the atoms is isotropic 
(x and z data overlapp) and the corresponding temperature is about 3.7 (2) /iK. In Fig. 
E^b) the cloud size as a function of time is given for an experiment in which the atoms 
were released from the crossed beam trap (off-set foci) and then allowed to evolve for 6 
ms inside the z beam before all potentials were switched off and the free expansion began. 
The asymptotic velocity for the elongated trap extracted from data given in Fig. Efb) gives 
T z = 0.43(2) /iK and T x = 4.4(2) /iK. The higher temperature in the x direction is caused by 
the focusing of the cloud and the extra collisions atoms undergo in the z direction because 
of the compressed cloud. The ratio is T x /T z = 10.2 which is quite close to the square of the 
mean aspect ratio = 3.3 for the collisionless expansion. Figure E^c) shows the ratio of the 
x to z size of the cloud at different times. Figure U\ shows the atomic cloud's image for 4 
different times during the free expansion. Figure [7^ a) shows the cloud at the crossed beam 
without expansion and Fig. Efb) is the same cloud after 4 ms free expansion. Fig. [TJc) 
is taken 6 ms after turning the vertical beam off and Fig. Efd) is the same cloud after 4 
ms of free expansion. The accurate simulation of the individual atomic trajectories for the 
free expansion is very time consuming because of the two-body cold collisions. In order to 
overcome this problem we emplo y th e direct simulation Monte Carlo (DSMC) method which 
was initially developed by Bird |26[ to simulate molecular gas dynamics. This method has 
been used for direct simulation of evaporative cooling (27 1 and free expansion of the cloud 



of atoms 



24j. Using this method the behavior of a cloud of atoms with 2 /iK temperature 



initially displaced 0.5zr from the potential minimum created by a gaussian beam with 25 
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FIG. 7: Images of the atomic cloud, (a) shows cloud at the beginning of the time of flight released 
from two crossed beams. In (b) the same cloud after 4 ms free expansion, (c) released from the 
crossed beam FORT 6 ms after turning off the x beam and (d) is the same cloud imaged 4 ms after 
the x beam has been extinguished. Each pixel in the images is 6 x 6 /im. 

fim beam waist was simulated. Figure El shows the simulated momentum distribution in 
the z (narrow distribution) and x (broad distribution) directions after 5 ms free expansion. 
The initial density was taken 1.8 x 10 14 atom cm~ 3 . As one can see the axial distribution 



deviates 
cloud [24 1 . 



ittle from a perfect gaussian. This has been reported previously for an expanding 



V. CONCLUSION 



In this paper we have studied the effect of geometrical changes of the trap on the atomic 
cloud evolution. We experimented with two configurations of crossed laser beams. In the 
first set of experiments, the foci of two orthogonal beams were positioned to conincide. The 
foci of the beams were subsequently offset for the second set of experiments. In each case, the 
behavior of the cloud upon abruptly switching off the x beam was studied. In the overlapped 
foci experiment we observed high atomic density domains, due to the strong anharmonic 
trapping potential. To reach the collisionless regime the trap density was reduced by lowering 
the power in the beams. We observed oscillations in the axial z direction with an amplitude 
up to a factor of 6 times the original cloud size. This implies that the momentum had been 
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FIG. 8: Simulated momentum distributions in z and x directions after 5 ms free expansion. Broken 
lines are the gaussian fit to the data. The narrow and broad curves are the z and x distributions 
respectively. 

reduced by a factor of 36. When the beams' foci were offset, cloud sizes as much as 10 
times the initial cloud size were achieved. In this later case we observed the atomic cloud 
to collapse in the z direction when released from certain positions in the trap. An energy 
transfer from the z to the x direction was observed for this release condition for clouds with 
high densities. 

Considering the fact that the cloud's center of mass gains a z velocity this system could 
:>e a useful tool in the study of quantum reflection of cold atoms from material surfaces 
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301 ] . Another possible avenue for further investigations would be subjecting a Bose- 
Einstein condensate to these potentials. 
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